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ABSTRACT 

Structural and optoelectronic properties of BaThOj cubic perovskite are calculated using all 
electrons full potential linearized augmented plane wave (FP-LAPW) method. Wide and direct 
band gap, 5.7 eV, of the compound predicts that it can be effectively used in UV based 
optoelectronic devices. Different characteristic peaks in the wide UV range emerges mainly due 
to the transition of electrons between valance band state O-p and conduction band states Ba-d, 
Ba-f, Th-fand Th-d. 
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1. Introduction 

Knowledge of physical properties of materials for their possible applications has always 
been a prime field of interest. Materials scientists even in today's technologically advanced 
society are still struggling for chief and efficient materials. Perovskites family contains a large 
number of compounds ranges from insulators to superconductors and from diamagnetic to 
colossal magneto-resistive (CMR) compounds. It is also expected that the ideal materials for 
spintronics applications could be double perovskites. 

Perovskites make huge portion of the mantle of the earth crust and therefore the 
investigation of physical properties of these compounds is highly desirable. They reveal many 
fascinating properties from both theoretical as well as experimental point of view. High 
thermoelectric power, ferroelectricity, superconductivity, charge ordering, spin dependent 
transport, colossal magneto-resistance, and the interplay of structural, magnetic and optical 
properties are commonly observed features in these materials. These materials are frequently 
used as sensors, substrates, catalytic electrodes in fuel cells and are also promising candidates for 
optoelectronics [1-9]. 

Barium thorate, BaThC>3, is one of the products of fission reactor and due to its 
importance in nuclear reactor Moreira et al. [10] studied its structural properties and Mishra et al. 
[11] studied its thermodynamic stability. Other physical properties of BaThC>3 like heat capacity 
[12] and Gibbs free energy [13,14] are also reported. Though barium thorate is an important 
member of the Perovskite family but even then no theoretical work has been published yet. 



In the present article, structural and optoelectronic properties of BaThC>3 cubic perovskite 
have been investigated. Lattice constant, bulk modulus, derivative of bulk modulus, band gap, 
density of states, complex dielectric function, complex refractive index, reflectivity, optical 
conductivity, energy loss function and absorption coefficient of the compound are calculated for 
the first time using self consistent full potential linearized augmented plane wave (FPLA-PW) 
method with Wu-Cohen generalized gradient approximation (GGA). It is expected that the 
present work will help in the understanding of optoelectronic behavior of BaThC>3 and the article 
will also cover the lack of theoretical data on the structural, electronic and optical properties of 
this interesting compound. 

2. Theory and computation 

Structural parameters of BaThOi are evaluated by fitting the unit cell energy versus unit 
cell volume using Birch-Murnaghan's equation of state [15]: 
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where E is the total energy of the unit cell, V is the unit cell volume, 5 and Bq are the bulk 
modulus at zero pressure and its derivative with pressure. Real and imaginary parts of dielectric 
function, ^(co) and £ 2 ( <u )> are calculated by the following relations [16, 17]: 
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where P ni i(k) is the dipole matrix element between initial and final states, S k is an energy 
surface with constant value, a) n? ^(A:) is energy difference between two states and p denotes the 



principal part of the integral. Real and imaginary parts of dielectric function are used to calculate 
refractive index and extension coefficient using relations: 

n(6>) = [KO) 2 + s 2 ^Y } 1 /2+e 1 (co) ]' 2 (4) 
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Real and imaginary parts of refractive index are further used in the following relation to evaluate 
normal incident reflectivity of BaThC^: 
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Similarly energy loss function L(o>)[18], absorption coefficient a(oj) and frequency dependent 
optical conductivity a(a>) are calculated by the following relations: 
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where VK ev is transition probability per unit time. 

In the present density functional calculations full potential linearized augmented plane 
wave (FP-LAPW) method with Wu-Cohen generalized gradient approximation [19], 
implemented in the wien2k code [20], is used to solve Kohan-Sham equations [21] for BaThC>3. 
In the full potential scheme wave function, potential and charge density are expanded into two 
different basis. The wave function is expanded in spherical harmonics in the atomic spheres 
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while outside the spheres (interstitial region) it is expanded in plane wave basis. The potential is 
also expanded in the same manner: 



where Eq. 10(a) is for inside and 10(b) is for outside of the atomic sphere. Inside the sphere the 
maximal value of / for the wave function expansion is l max =10 and is spherically symmetric 
while outside the sphere it is constant. R M t is chosen in such a way that there is no charge 
leakage from the core and the total energy convergence is ensured. R M t values of 2.5, 2.19, 1.94 
a.u. are used for Ba, Th, O respectively. For wave function in the interstitial region the plane 
wave cut-off value of K max = II R M t is chosen. 2000 k points are used in the Brillouin zone 
integration and convergence is checked through self consistency. The perovskite structure of 
BaThC>3 is shown in Fig. 1 . Ba atom in the unit cell is at (0,0,0), Th atom at (0.5,0.5,0.5) while O 
atoms at the faces of the unit cell, (0.5,0.5,0), (0.5,0,0.5), (0,0.5,0.5). 
3. Results and Discussion 

Structural and optical properties of BaThC>3 are calculated by all electrons ab-initio full 
potential linearized augmented plane wave (FP-LAPW) method. The unit cell volume of BaThC>3 
is varied in the neighborhood of the experimental value 4.48 A [10] and the corresponding 
energies are calculated using Wu-Cohen GGA scheme. Birch-Murnaghan's equation of state [15] 
is used to plot, relation between energy and volume shown in the Fig. 2. From the plot relaxed 
lattice constant is obtained by taking the volume corresponding to the minimum energy. The 
optimized parameters are quoted in Tableland are compared with the other available data. It is 
evident from the table that our calculated lattice parameter is in good agreement to the 
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experimental [10] and other theoretical results [31, 32]. The large value (124 GPa) of bulk 
modulus suggest that BaThC>3 is less compressible with pressure derivative of 2.41 . 

Band structure and density of states are calculated with the well converged self consistent 
solution of FP-LAPW, shown in Figs. 3 and 4. It is clear from Fig. 3 that, BaThC>3 is a direct 
band gap compound with a band gap of 5.7 eV at the T symmetry point. The calculated direct 
band gaps at other symmetry points are also presented in Table 1 . Wide and direct band gap at 
the T symmetry point has also been reported for similar compounds like BaTiC>3 [22], BaZrC>3 
[23] and BaHfC>3 [24] which confirms the correctness of our band structure calculation for 
BaTh0 3 . 

The calculated density of states for BaThC>3 is shown in Fig. 4. It is clear from the figure 
that the density of states can be divided into three regions. The first region is from -13.57 eV to 
-1 1.48 eV with a major contribution from Ba-5p states while small contribution from Th-6p 
states. The next region is valence band, ranges from -3.7 eV to 0.0 eV and the main contribution 
is due to 0-2p states. The last region is conduction band, which ranges from 3.33 eV to 13.66 
eV. In this band, the prominent contribution comes from Ba-4f states with a noticeable part from 
Ba-5d, Th-5f and Th-6d states. 

In order to envision the bonding nature and charge transfer the charge density in the (100) 
and (110) planes is plotted in Fig. 5(a, b). From Fig. 5(a), it can be seen that there is a large 
transfer of charge among Ba and O atoms with a very small contour extension among the ions 
and hence Ba-0 bond is strongly ionic with very weak covalent nature. Fig. 5(b) shows that 
electrons are strongly shared and distributed along the Th-0 bond. Hence, their bond is strongly 
covalent with a bond length 2.27 A . Similar charge distribution and bonding nature has been 
reported for a similar compound, BaHfC>3 [24] . 
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Optical spectra of a compound provide useful information about the internal structure of 
that material. Keeping in view the importance of optical properties of BaThC>3, different 
frequency dependent optical parameters are described in terms of electronic properties to 
recognize the optoelectronic nature of the compound for different devices applications. 

The dielectric function,£(a>) = ^(co) + ie 2 (joS), describes the optical response of the 
medium to the incident photons, E = hco. The imaginary part of the dielectric function, £ 2 (w), is 
the most important parameter of the optical properties of a material and therefore its details study 
is important. The variation of the imaginary part of the dielectric function with the incident 
photon energy is shown in Fig. 6. It is clear from the figure that the critical point of £ 2 ( a) ) is 
around 5.3 eV. Two high peaks A and C at 8.7 eV and 14.9 eV along with some noticeable peaks 
B, D, E, F and G are located at 12.5, 17.9, 18.5, 20 and 26 eV respectively. The origin of these 
peaks lies in the inter band transitions which can be related to the density of states of the 
compound shown in Fig. 4. It can be seen in the figure that the peaks A, B, C and D are due to 
the transition of electrons from O-p states to Ba-d, Th-f, Ba-f and Th-d states respectively. While 
the peaks E, F and G emerges due to the transition of electrons from mixed Ba-p and Th-p states 
below the valence band to Th-F, Ba-f and Th-d states in conduction band respectively. 

The optical band gap obtained from the imaginary part of the dielectric function is 5.6 eV 
which is close to the band gap obtained from the band structure (5.7 eV). It is well understood 
that materials with band gaps larger than 3.1 eV works well in the ultraviolet (UV) region of the 
spectrum [25-29]. Hence this direct and wide band gap material could be suitable for the high 
frequency UV devices applications. As the compound also shows strong absorbing nature in the 
different parts of spectra, ranges 7- 27 eV. So, it can also be used as a filter for various energies 
in the UV spectrum. 



Different important frequency dependent optical parameters like; real part of dielectric 
function £i(co), extinction coefficient /c(o>), absorption coefficient a(o>), energy loss 
function L(a>), refractive index n(a>), reflectivity R (o>) and optical conductivity cr(a>)are plotted 
in Fig. 7. The plot of the real part of the dielectric function (Fig. 7a) explains different important 
physical aspects of BaThC>3. It is clear from the figure that zero frequency limit, £i(0), is 2.95. 
The magnitude of £i(o>) increases from zero frequency limit to peak value 7.9 at 6.9 eV. After 
maximum value it starts decreasing and becomes flatten with small variations. It is further noted 
that the real part of dielectric function becomes negative in the energy ranges 15.19-15.6 3 eV, 
16.01- 16.39 eV, 18.65- 19.06 eV, 20.01- 23.68 eV and 26.62-28.15 eV. As materials behave 
metallic for negative values of £i(o>) and are dielectric otherwise [30]. Therefore for these 
ranges of energy BaThC>3 is metallic. Interestingly, the refractive index (Fig. 7e) follows the 
trend of £i(o>). The zero frequency limit of the refractive index is 1.7 with a maximum value of 
2.7 at 8.1 eV. 

The figure also reveals similarities in the trends of extinction coefficient k(ui) and 
absorption coefficient a(cS) (Fig. 7b and 7c). The Gaussian smearing value in the present 
calculations is 0.1 eV. The critical point of /c(o>) is 6.3 eV and of <x(oj) is 6.8 eV. Both have 
strong response to the incident photons in the range 8-29 eV. The highest peak value of 
extinction coefficient is 8.7 and of absorption coefficient is 218 x 10~ 4 cm _1 at 20.2 eV. 

Frequency corresponding to the plasma resonance can be calculated from the energy loss 
spectra, shown in Fig. 7(d). The maximum resonant energy loss is at 28.5 eV which corresponds 
to plasma frequency 2.7 x 10 15 Hz. From Fig. 7(f), the zero frequency limit of reflectivity for 
BaThC>3 is found to be 0.065. There are high reflection peaks at energies 8.7 eV, 15.25 eV, 22.1 
eV and 27.8 eV corresponding to the negative values of £i(o>) . The optical conductivity (Fig.7g) 
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starts responding to the applied energy field from 6.5 eV. The good response is found in the 
range 7 to 28eV. Maximum optical conductivity of the compound is at 14.87 eV of 
magnitude 9288/2 -1 cm' 1 . 

4. Conclusion 

Density functional FP-LAPW calculations are performed on BaThC>3 for the first time to 
investigate structural and optoelectronic properties. The calculated lattice constant 4.55 A is in 
good agreement with the available data. The calculated band structure predicts that the 
compound has a wide and direct band gap of 5.7 eV. On the basis of the direct wide band gap 
and the spectra of imaginary part of the dielectric function it can be concluded that this material 
will be useful for optoelectronic devices in the UV region of the spectrum. 
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Figure Captions 
Fig.l: Crystal structure of cubic perovskite BaTh03. 
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Fig.2: Variation of total energy as a function of unit cell volume of BaThCh 
Fig.3: Band structure of BaThC>3. 
Fig.4: Total density of states of BaThC>3. 

Fig.5: Total electron density for BaThC>3: (a) in the (100) plane (b) in the (110) plane 
Fig.6: Imaginary part of dielectric function as a function of energy for BaThC>3. 
Fig.7: Frequency dependent optical parameters for BaTh03. 
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Table 1: Lattice constant, bulk modulus, pressure derivative of bulk modulus, Band gap of 
BaTh0 3 . 



This Work Experimental work Other work 



Lattice constant o(A) 4.55 4.48 1 4.43 11 , 4.53 m 

Bulk modulus B(GPa) 124.34 

Derivative of bulk modulus B 2.41 

Band gap £^ r (eV) 5.70 

Band gap £"£~ R (eV) 7.70 

Band gap E™~ m (qV) 7.90 

Band gap ff%~ x (eV) 7_A0 
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